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Inside the nucleus, DNA replication is organized at 
discrete sites called replication factories, consisting of 
DNA polymerases and other replication proteins. Rep- 
lication factories play important roles in coordinating rep- 
lication and in responding to replication stress. However, 
it remains unknown how replicons are organized for pro- 
cessing at each replication factory. Here we address this 
question using budding yeast. We analyze how individ- 
ual replicons dynamically organized a replication factory 
using live-cell imaging and investigate how replication 
factories were structured using super-resolution micros- 
copy. Surprisingly, we show that the grouping of replicons 



within factories is highly variable from cell to cell. Once 
associated, however, replicons stay together relatively 
stably to maintain replication factories. We derive a coher- 
ent genome-wide mathematical model showing how neigh- 
boring replicons became associated stochastically to form 
replication factories, which was validated by independent 
microscopy-based analyses. This study not only reveals 
the fundamental principles promoting replication factory 
organization in budding yeast, but also provides insight 
into general mechanisms by which chromosomes orga- 
nize sub-nuclear structures. 



Introduction 

DNA replication initiates at replication origins and proceeds 
as replication forks move along parental DNA. At each replica- 
tion fork, DNA polymerases and accessory proteins, such as 
proliferating cell nuclear antigen (PCNA), form a large com- 
plex (called the replisome) that conducts de novo DNA synthe- 
sis. Current evidence suggests that in eukaryotes, the two sister 
replisomes initiated from each origin stay associated (Falaschi, 
2000; Kitamura et al., 2006; Ligasova et al., 2009) and normally 
coordinate replication of the entire replicon (Conti et al., 2007; 
Natsume and Tanaka, 2010). Nevertheless, if one fork stalls, for 
example due to DNA damage, the other fork can still continue 
DNA synthesis (Doksani et al., 2009; Yardimci et al, 2010). 
In eukaryotes, multiple pairs of sister replisomes are grouped 
together at discrete sites called replication factories (Berezney 
et al., 2000; Kitamura et al, 2006; Gillespie and Blow, 2010). 
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Upon replication stress or DNA damage, a replication factory 
defines an important boundary, inside of which dormant origins 
can initiate replication and thus evade checkpoint-dependent 
inhibition (Dimitrova and Gilbert, 2000; Ge and Blow, 2010; 
Thomson et al., 2010). Moreover, replication factories may facili- 
tate efficient DNA synthesis by creating high concentrations of 
replication proteins (Natsume and Tanaka, 2010). 

However, in spite of their functional importance, principles 
that underpin the organization of replication factories remain 
elusive. There has been a long-standing debate about whether 
replication factories are organized by deterministic constraints 
or by stochastic processes. Several reports have suggested that 
the organization is facilitated by particular regulatory factors 
(Yan et al., 1998) or post-translational modifications (Rossi et al., 
1999). Nevertheless, these observations still leave open whether 
or not replicons are grouped together stochastically. 
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Results and discussion 

Live-cell microscopy reveals that: 
neighboring replicons often form replication 
factories but their grouping is highly 
variable from cell to cell 

To investigate the organization of replication factories, we ana- 
lyzed the behavior of replicons using live-cell microscopy in the 
budding yeast Saccharomyces cerevisiae. We chose a region on 
chromosome VII with four adjacent replicons (Fig. 1 A). We 
selected one locus on each replicon such that all four loci show 
the same average replication timing (Yabuki et al., 2002). We 
integrated tetO and lacO arrays at paired loci in three different 
yeast strains (strains #1-3; Fig. 1 A, bottom). These arrays 
bound TetR and LacI proteins, fused with cyan and green fluor- 
escent proteins (CFP, GFP), respectively, and were thus visual- 
ized as small fluorescent dots. The fluorescent dots increased in 
intensity upon DNA replication as the number of arrays was 
doubled, thus defining their replication timing by microscopy 
(Kitamura et al., 2006). Fluorescent dots did not increase their 
intensity if DNA replication was inhibited. 

It is known that there is cell-to-cell variation in the time 
that any particular origin initiates (Bechhoefer and Rhind, 
2012), so that the two marked loci in each strain did not always 
replicate with similar timing (Fig. SI A). To analyze how repli- 
cons are gathered into factories, we focused on the cells in which 
the two marked loci replicated with similar timing (difference 
<3 min). We then determined whether the two loci replicated in 
close spatial proximity (<350 nm apart for >2 min), meaning 
that they replicated in the same factory (Fig. 1, B [top] and C). 
Alternatively, they could replicate in different factories (Fig. 1, 
B [bottom] and D). The same proximity threshold was previously 
used to define the association of sister replisomes (Kitamura et al., 
2006) and was independently validated as described below. 

Using this protocol, the two marked loci in strain #1 that 
are 70 kb apart (the distance from mid-tetO to mid-/ac0) repli- 
cated in the same factory in 43% of cells (10/23) and in different 
factories in 57% of cells (13/23; Fig. 2 A), suggesting that 
grouping of replicons within factories can vary from cell to cell. 
In contrast, in strains #2 (loci 131 kb apart) and #3 (loci 192 kb 
apart), the two marked loci replicated in the same factory less 
frequently (Fig. 2 A and Fig. SI, B-D): 11% (2/19) and 5% (1/19), 
respectively (P = 0.037 and 0.0059, respectively, compared with 
strain #1). Thus, replicons close along a chromosome were 
often processed for replication in the same factory, but replicons 
farther apart replicated more frequently in different factories. 

We next addressed if the two marked loci being close to- 
gether before replication leads to a higher probability of them 
replicating in the same factory. We determined the distance be- 
tween the loci in strain #1 before replication, and showed that its 
distribution had no correlation with whether the loci were repli- 
cated in the same or different factories (Fig. 2 B and Fig. S2 A). 
In contrast, when we measured the distance between loci soon 
after replication (until +11 min when the fork that passed 
through the tetO array is expected to have terminated; Fig. 1 A), 
loci that replicated in the same factory tended to remain closer 
than loci that replicated in different factories (Fig. 2 C; refer to 



a control in Fig. S2 B). If the same measurement is made >1 1 min 
after replication, the difference in the distance between loci 
replicating in the same and different factories was again lost 
(Fig. 2 D). This suggests that the factory configuration per- 
sists for a few minutes after loci have replicated, but is lost 
within 1 1 min (perhaps when replication of a replicon is com- 
pleted). This is consistent with a relatively high rate of diffu- 
sion of chromosomal loci within the nucleus (see Materials and 
methods, Mathematical modeling ii). Moreover, the difference 
between the same versus different factory groups (Fig. 2 C) 
validates the original loci-proximity threshold to estimate repli- 
cation in the same factory, as an inaccurate threshold would not 
create such a difference. 

Super-resolution microscopy reveals the 
distribution of replicons within individual 
replication factories 

Using live-cell imaging, we previously visualized replication fac- 
tories as bright globular foci of GFP-PCNA, which showed 
dynamic changes in shape and location during S phase (Kitamura 
et al., 2006). To provide a better understanding of replication fac- 
tory organization, we observed GFP-PCNA by super-resolution 
structured-illumination microscopy (Schermelleh et al., 2008). 
Cells were fixed and analyzed at different times after release 
from a-factor treatment (time 0), allowing us to monitor factory 
number and fluorescence intensity at different times. This showed 
that the number of replication factories increased to a peak value 
of 73 ± 8 (mean ± SD) in mid-S phase (Fig. 3, A and B). 

We next evaluated the number of replisome pairs present in 
each replication factory. Published replication profiles showing 
the replication timing of the whole genome (Raghuraman et al., 
2001 ; Yabuki et al., 2002) are an average from a large number of 
cells and do not accurately represent replication in individual cells. 
For example, the number of "origin" peaks found in an averaged 
profile is usually larger than the number of replication initia- 
tion events in individual cells because replication origins initiate 
with <100% efficiency (de Moura et al., 2010). To estimate the 
total number of forks, we simulated the replication in individual 
cells, using a technique originally established for chromosome VI 
(de Moura et al., 2010) and subsequently extended to the whole 
genome (unpublished data), excluding the ribosomal DNA (rDNA) 
region. Using this simulation, we estimated that 242 ± 24 (mean ± 
SD) replication forks were present at the peak of DNA replication 
(Fig. 3 C). Using published data (Linskens and Huberman, 1988; 
Pasero et al., 2002), we estimated that 60 replication forks were 
present, on average, in the rDNA region. We estimated that on 
average, 302 forks (242 + 60) are found in the whole nucleus at the 
peak of DNA replication. We then assigned these 302 replication 
forks to each of the replication factories, assuming that the inte- 
grated GFP-PCNA signal in each factory is proportional to the 
number of forks it contains (Fig. 3 D). In this way, we were able to 
estimate the number of replisome pairs (at sister replication forks) 
present in each replication factory (Fig. 3 E). Intriguingly, our result 
suggests that many replication factories (39%) consist of just one 
pair of sister replisomes, as was previously predicted (Berezney 
et al., 2000). Factories with more replisome pairs are less abun- 
dant, with only 17% of factories having >4 replisome pairs. 
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Figure 1 . Live-cell imaging reveals assembly of replicons at replication factories. (A) Three yeast strains were constructed to analyze assembly of three 
different replicon pairs at an individual factory. In strains #1 (T7278], #2 (T9184), and #3 [17277] with TetR-3xCFP and GFP-Lacl, the indicated two 
chromosomal loci were marked by integration of fe/Ox224 and /acOx256. The replication timing profile (top) was obtained from Yabuki et al. (2002); 
0 min represents the time of release from a-factor arrest. (B) Models for replication of two replicons at the same factory and at different factories. (C and D) 
Representative examples of strain #1 showing replication of CFP and GFP dots at the same factory (C) and at different factories (D). Cells were released 
from a-factor treatment. CFP and GFP images were acquired every 1 min for 65 min. The intensity of CFP/GFP dots and their distance were measured in 
23 cells of strain #1 in single experiments. The change in intensity of each dot was fitted by a sigmoidal dose-response curve (see Materials and methods; 
see R 2 , representing fitness, below) in which a midpoint in the increase of its intensity was defined as its replication time. Then 0 min was defined as the 
mid-replication time of two dots. R 2 = 0.70 (C, CFP), 0.81 (C, GFP), 0.78 (D, CFP), and 0.74 (D, GFP). Bars, 1 urn. 
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Figure 2. Assembly of replicons for a replication factory varies from cell to cell and closer ones show more frequent assembly. (A) In each strain of # 1 -3 
(see Fig. 1 A), 1 9-23 cells were analyzed as in Fig. 1 , C and D, and categorized based on their two fluorescent dots being replicated at the same (yellow) 
or different factories (light green). 0 min was defined as in Fig. 1 , C and D. Replication at the same factory was evaluated by close proximity (<350 nm 
for 2 min or longer) of the dots within the time window from —3 min to +2 min (orange rectangle). (B-D) The distance between the two dots was measured 
at each time point in 23 cells of strain #1 , as shown in A. Its distribution was then plotted as cumulative frequency, separately for replication at the same 
(red line) and different factories (blue line), in three time windows; i.e., (B) before replication, (C) soon after replication (excluding the time window —3 min 
to +2 min for the evaluation; see A), and (D) > 1 1 min after replication, n (number of time points) = 1 90 (B, same), 254 (B, different), 87 (C, same), 1 1 7 
(C, different), 140 (D, same), 182 (D, different). 
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Figure 3. Super-resolution microscopy reveals the organization of replication factories. (A and B) Replication factories observed by super-resolution micros- 
copy. Cells (T8375) with GFP-POL30 (PCNA), SPC42-mCherry (a component of spindle pole body, SPB), and NIC96-mCherry (a component of the nuclear 
pore complex, NPC) were released from a-factor treatment (defined as 0 min). (A) A bright-field image, a fluorescence image (GFP, green; mCherry, red), and 
its 3D rendering in a representative cell, which was fixed at 40 min. Yellow arrows on the 3D rendering image indicate representative replication factories that 
are estimated to contain a single replisome pair (see D). Bars: (bright-field image) 2 urn; (fluorescence image) 0.5 urn. (B) The number of PCNA foci (mean ± 
SD) within the nucleus along the time course. Red line indicates the peak of DNA replication. (C) The number of replication forks along the time course (0 min; 
DNA replication initiated in 50% cells), based on the in silico analysis. The blue dots represent the mean fork numbers with 1-min intervals (green broken lines: 
± SD). The number of forks reached the maximum (242 ± 24, mean ± SD) at 30 min (red line). The data did not include the forks along the rDNA region. 
(D) Schematic diagrams explaining how the number of replisomes (at replication forks) was estimated at each factory. This number was binned to even integers 
as sister replisomes associate (right; Kitamura et al., 2006). (E) Relative fractions of replication factories containing each number of sister replisome pairs. The 
intensities of individual replication factories were analyzed at the peak of replication (40 min in B) as explained in D, collectively in 23 cells. (F) Representative 
super-resolution fluorescence images, with schematic diagrams of interpretation. Cells (T8659) with TagBFP-POL30 (PCNA) and ORC2-mCherry were treated 
as in A. At 45 min, cells were labeled with the thymidine analogue EdU for 12 min, followed by fixation. Bars: (left) 0.5 pm; (right) 0.15 urn. 
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If the majority of replication factories comprise one or 
two pairs of sister replisomes, as shown in Fig. 3 E, we may be 
able to visualize their association with nascent DNA using 
super-resolution microscopy. To observe both factories and rep- 
lication origins, we fused PCNA and origin recognition complex 
subunit 2 (Orc2) with different fluorescent proteins (Fig. 3 F; 
blue and red, respectively). We also visualized nascent DNA 
by a 12-min pulse of 5-ethynyl-2'-deoxyuridine (EdU, green). 
To simplify our interpretation, we focused on factories that were 
well separated from others and were closely associated with 
Orc2 and EdU. We observed a high proportion of cases (22/40; 
P = 0.0037) where the labels were present in the order PCNA- 
EdU-Orc2 (Fig. 3 F, top). This may represent a replicon emerg- 
ing from a factory during replication, with nascent DNA lying 
behind the factory, and the origin (marked by Orc2) behind the 
nascent DNA. In addition, we found several examples that were 
consistent with two replicons coming out of a single factory 
(Fig. 3 F, bottom). We also studied the direction in which repli- 
cated DNA came out of a factory, and found it to be random rel- 
ative to the factory-nucleus center axis (Fig. S3 A). 

Genome-wide stochastic assembly of 
replicons recapitulates the organization 
of replication factories within the 
whole nucleus 

How are replisome pairs assembled into replication factories? 
As shown in Fig. 2 A, this assembly differed from cell to cell. 
Given this, the simplest hypothesis is that factory assembly is 
stochastic. If replisome pairs are dispersed in the nucleus and 
they randomly assemble into a factory, the number of replisome 
pairs per factory should follow a Poisson distribution (Motulsky, 
2010). Consistent with this, our estimate of the number of repli- 
somes in each factory (Fig. 3 E) approximately matched a Poisson 
distribution (Fig. 4 A). 

Prompted by this finding, we developed a simple mathe- 
matical model for factory formation (Fig. 4 B; see Materials and 
methods, Mathematical modeling i and ii) where two adjacent 
pairs of sister replisomes, connected by a chromosome region 
(length d), are represented by two particles on a string. The two 
particles diffusing rapidly within the nucleus can aggregate with 
a binding energy J if they come within distance e from each 
other. Using basic statistical physics, we derived the probability 
P{d) of the two particles being associated as: 

W Ad 3 +] 

where A is a constant depending on J, s, and temperature. We 
derived the parameter A as 8.7 x 10~ 6 kb~ 3 by fitting this equa- 
tion (Fig. 4 C and Fig. S3 B) to the results of replicon assembly 
frequency, obtained in Fig. 2 A. The model fitted the data well 
(R 2 = 0.99; Fig. 4 C). If the diameter of a sister replisome pair 
is 90 nm (i.e., e = 90 nm; Baddeley et al., 2010), / is about 
— 12.5 kj/mol (—5.1 k B T) at ambient temperature, which is an 
energy involved in a typical weak protein-protein interaction 
(Baxter et al., 1998; Rippe, 2007). 



Next, we extended this model to the whole genome by 
combining it with the mathematical simulation of the replica- 
tion profiles in individual cells (de Moura et al., 2010; unpub- 
lished data). Note that all parameters in this simulation were 
determined from the known profiles of DNA replication and we 
did not tune any of the parameters. For each cell in one million 
simulations we acquired a "snapshot" of the replication fork 
positions on the genome at the peak of DNA replication (Fig. 4 D, 
top; see Materials and methods, Mathematical modeling iii). 
The collection of these snapshots gave the distribution of the 
distance between neighboring sister fork pairs, as shown in 
Fig. 4 E. We then used our equation for P{d) as described above 
(Fig. 4 C) to determine the probability of fork pairs aggregat- 
ing into a factory given their distance d along a chromosome 
(Fig. 4 D, bottom). This gave an "in silico" distribution of num- 
bers of sister fork pairs per replication factory (Fig. 4 F, green 
squares). Remarkably, the in silico result was very similar to the 
microscopy observations (Fig. 4 F, red circles). Thus, from the 
frequency that adjacent sister replisome pairs associate with 
one another (Fig. 2 A) we were able to accurately recapitulate 
the genome-wide replicon distribution in replication factories 
(Fig. 3 E) by assuming stochastic assembly of replicons. 

Implications of this study 

Our results suggest that it is mainly neighboring replicons on 
a chromosome that are brought together in factories (Figs. 2 A 
and 4 C). This conclusion is widely anticipated (Jackson and 
Pombo, 1998; Berezney et al., 2000; Gillespie and Blow, 2010) 
and is consistent with the observed clustering of active replicons 
on DNA fibers (Tuduri et al., 2010) and the high rate of associa- 
tion of neighboring DNA sequences observed in chromosome 
conformation capture assays (Duan et al., 2010). 

Previous data suggested that sister replisomes stay associ- 
ated with each other during replication (Kitamura et al., 2006). 
Do the sister replisome pairs present in the same factory remain 
continuously associated with each other, or do they often disso- 
ciate? We found that the physical distance between two marked 
loci tends to remain shorter after they were replicated in the 
same factory than after being replicated in different factories 
(Fig. 2 C). This suggests that replicons are stably associated for 
a significant period once they are brought together. Furthermore, 
the association between replicons fitted very well a stochastic 
thermodynamic process (Fig. 4, C and F) in which their associa- 
tion represents a low energy state. Thus, the association between 
replicons is expected to be relatively stable. 

Our results suggest that a replication factory is mainly 
organized by associated neighboring replicons on the same 
chromosome. However, there may be other types of replicon 
association that contributes to factory formation, although they 
are unlikely to constitute a major population. For example, as 
centromeres cluster at a spindle pole in budding yeast (Duan 
et al., 2010), replicons at centromeres on different chromosomes 
may be associated to organize a factory (Natsume et al., 2013). 
Moreover, a DNA binding protein that recognizes a common 
DNA motif may bring together a subgroup of replicons that are 
distant on the same chromosome or present on different chro- 
mosomes (Knott et al., 2012). 
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Figure 4. Stochastic assembly of sister replisome pairs recapitulates their observed distribution at individual replication factories. (A) Distribution of the 
number of sister replisome pairs per replication factory (Fig. 3 E) approximately fits a Poisson distribution (X = 1 .60, R 2 = 0.98). (B) Schematic diagram of 
"Particles on a string" model. See detail in Materials and methods, Mathematical modeling i. (C) We fitted the parameter A in P(d) = l/(Aa° + 1), using 
the replicon assembly frequency result, obtained in strain #1 and #2 (Fig. 2 A; 23 and 19 cells were analyzed, respectively). Fig. S3 B explains how the 
chromosomal distance between sister replisome pairs was obtained in strain #1, #2, and #3. (D) Schematic diagram showing the example of positions 
of sister replisome pairs along a chromosome (top), and their stochastic assembly (bottom) based on the grouping probability shown in C. (E) Distribution 
of the distance (replicated DNA is not counted, as it is looped out; see C) between neighboring sister replisome pairs along a chromosome, obtained 
from the simulation. Relative fractions of the pairs at the indicated distance (each 5-kb window) were obtained from one million simulations, at the peak 
of DNA replication (at 30 min in Fig. 3 C). Fractions for >250 kb are shown together at 250 kb. (F) Distributions of the number of sister replisome pairs 
per replication factory: comparison between in vivo data (obtained in Fig. 3 E) and in silico data (as shown in B-E). Pearson correlation coefficient 
r= 0.9953. P < 0.0001. 
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It is intriguing to consider whether replication factories are 
organized in mammalian cells in similar or different ways com- 
pared with budding yeast. Some mammalian factories formed in 
late S phase are as large as the whole yeast nucleus (see Fig. 3 in 
Natsume and Tanaka, 2010). A chromosome region, replicated 
at such a large mammalian factory, represents a chromosome 
territory — a stable unit maintained until the next cycle (Jackson 
and Pombo, 1998). It is unlikely that a chromosome region repli- 
cated at a yeast factory represents such a stable unit because the 
motion of a chromosome locus (on a chromosome arm) covers a 
large proportion of the yeast nuclear volume (see Materials and 
methods, Mathematical modeling ii). Nonetheless, a mammalian 
chromosome territory is composed of smaller units called topo- 
logically associating domains (TADs), inside of which chromo- 
somes show more dynamic behaviors (Gibcus and Dekker, 2013). 
Stochastic assembly of yeast replicons forming replication facto- 
ries may correspond to the way that chromosomal DNA associ- 
ates within a mammalian TAD. Indeed, mammalian cells also 
have small factories, some of which contain only 1-3 replicons 
(Berezney et al, 2000), as we found in yeast. Moreover, the large 
factories seen in mammalian cells may actually be composed of 
several small ones (Leonhardt et al, 2000; Baddeley et al., 2010). 
Such small factories of mammalian cells may be organized in 
the same way as we found in yeast cells. 

A stochastic assembly mechanism may provide robustness 
to factory organization. It is relatively easy to establish — all that 
is required is that some replisome components have an affinity 
for another replisome component. In a deterministic assembly 
scheme, failure to incorporate one component might cause failure 
of the entire factory network, whereas in a stochastic scheme, 
each individual interaction is independent of the status of the 
others. This could be important in vertebrate cells in responding 
to replication stress when a replication factory defines the 
boundary, inside of which dormant origins can initiate and com- 
plete replication for the region between two stalled replication 
forks (Ge and Blow, 2010; Thomson et al, 2010). 

In addition to organizing DNA replication, replication fac- 
tories (foci) are likely to represent a fundamental feature of chro- 
mosome organization (Jackson and Pombo, 1998; Berezney et al., 
2000). Using budding yeast as a model organism, we have found 
that the replicons making up individual replication factories are 
highly variable from cell to cell. Our results suggest that neigh- 
boring replicons are assembled stochastically and stay associated 
together to maintain replication factories relatively stably. Our 
study gives an important insight not only to the organization of 
DNA replication within the nucleus, but also to general mecha- 
nisms by which chromosomes organize sub-nuclear structures 
such as transcription factories and repair foci (Lisby and Rothstein, 
2004; Sutherland and Bickmore, 2009). 

Materials and methods 

Yeast genetics and molecular biology 

All S. cerevisiae strains used in this study were MATa haploid strains 
with the W303 genetic background (trpl-1, Ieu2-3,1 12, ura3-l, his3- 
1 1 , 15, ade2-l, can 1-100). The methods in yeast genetics, a-factor treat- 
ment for cell synchronization (Amberg et al., 2005; Kitamura et al., 2006), 
the TetR-GFP/fef and GFP-Lacl//ac operator system (Straight et al., 1996; 



Michaelis et al., 1997), and the construct of TetR-3xCFP (Bressan et al., 
2004) were as described previously. Cells were cultured in YPA medium 
containing 2% glucose at 25°C unless otherwise stated. 

tetOx224 (1 1 .2 kb) and lacOx256 (10.1 kb) were integrated by a 
two-step "pop-in and pop-out" method (Struhl, 1983), at the chromosomal 
loci shown in Fig. 1 A, as follows. tetOx224 was integrated to chromo- 
some VII at 660847 (base pairs from the left telomere) within the replicon 
including autonomously replicating sequence, ARS727 (the ARS727 repli- 
con), in strains #1 , #2, and #3. lacOx256 was integrated to chromosome 
VII at 720007 within the ARS728 replicon, at 78 1 554 within the ARS729 
replicon and at 842709 within the ARS731 replicon in strains #1 , #2, and #3, 
respectively. To make constructs for integration, DNA fragment spanning 
~500 bp upstream (toward the left telomere) and ~1 kb downstream 
(toward the right telomere) from integration sites were amplified by PCR 
and cloned into the pRS306 plasmids containing URA3 as an auxotroph 
marker. The operator arrays were then inserted between these upstream 
and downstream genomic DNA fragments. The plasmids were subsequently 
cut within the upstream fragment and inserted into each locus. The cells 
with the plasmid backbone popped out were selected on 5-fluoroorotic 
acid-containing media (Amberg et al., 2005). Plasmid integrations and 
subsequent loss of plasmid backbones were confirmed by PCR amplifica- 
tions of relevant regions. 

When the N terminus of PCNA [POL30] was tagged with a fluores- 
cent protein at its original locus, the growth of haploid cells was severely 
retarded (Kitamura et al., 2006). We therefore integrated a single copy of 
PCNA, tagged at its N terminus, into an auxotroph marker TRP while 
PCNA is intact at its original locus. It has been shown in fission yeast that 
fluorescently tagged PCNA behaves similarly to untagged PCNA (Meister 
et al., 2007). PCNA and its promoter were first cloned into pRS404 
(pT1046) and then DNA fragments of yEGFP (from pKT128; Sheff and 
Thorn, 2004) and TagBFP (from pTagBFP-C; Evrogen) were inserted be- 
tween the promoter and the coding region to construct pT1056 and 
pTl 448, respectively. SPC42 and N1C96 were tagged at their C termini at 
the original gene loci with a one-step PCR method (Knop et al., 1 999), 
using 4xmCherry cassette from pT909 as a PCR template. The plasmid 
pT909 was constructed by multiplying the mCherry gene in pKS391 (Snaith 
et al., 2005). Strains with the tagged genes grew normally at temperatures 
used in this study. 

To facilitate EdU incorporation (Salic and Mitchison, 2008), five 
copies of the herpes simplex thymidine kinase gene were expressed from 
GPD1 promoters (Dahmann et al., 1995), and the human equilibrative 
nucleoside transporter 1 [ENT1] gene (with codon usage optimized for yeast) 
was expressed from ADH1 promoter (a gift from K. Shirahige, University of 
Tokyo, Tokyo, Japan). 

Live-cell imaging 

The procedures for time-lapse fluorescence microscopy were as described 
previously (Kitamura et al., 2006; Tanaka et al., 2010) unless otherwise 
stated. During image acquisition, cells were suspended in synthetic complete/ 
YPA (3:1 ratio) medium containing 2% glucose at 25°C (ambient tempera- 
ture). For image acquisition, we used a microscope (DeltaVision RT; Applied 
Precision), UPlanSApo 1 OOx/NA 1 .40 objective lens (Olympus), a CCD 
camera (CoolSnap HQ; Photometries), and SoftWoRx software (Applied 
Precision). CFP and GFP signals were discriminated with the 89006 ET 
filter set (Chroma Technology Corp.). We acquired nine z-sections (0.7 pm 
apart), which were subsequently analyzed with SoftWoRx, Volocity (Perkin- 
Elmer), and Image-Pro Plus (Media Cybernetics) software. For two-dimensional 
presentation in figures, z-sections were deconvoluted with SoftWoRx and 
projected to two-dimensional images. For quantitative analyses, images 
before deconvolution were used. 

Super-resolution structured illumination microscopy 

Cells were fixed with 2% paraformaldehyde for 30 min and washed with 
PBS. When cells were pulse-labeled with 1 mM EdU (for 12 min; Fig. 3 F), 
they were further processed for the Click-iT EdU reaction as described in 
the protocol CI 0337 (Invitrogen), followed by a PBS wash. The images 
were acquired with structured illumination microscopy (Schermelleh et al., 
2008) as described previously (Hattersley et al., 201 1). In brief, we used 
the OMX system (version 3; Applied Precision) equipped with 405-, 488-, and 
593-nm solid-state lasers. To acquire images, we used a 512 x 512 electron- 
multiplying charge-coupled device camera (Cascade II; Photometries) and 
a UPlanSApochromat 1 00x/l .4 NA lens. A coherent scrambled laser 
through a diffraction grating was used to generate the structured illumina- 
tion. We used SoftWorx to process raw images for reconstruction to reveal 
structures in three dimensions (Gustafsson et al., 2008). 
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image and data analyses 

The grouping of replicons into the same factory or different factories was 
evaluated based on the distance between two dots (CFP and GFP signals] 
during replication. For an accurate evaluation, we analyzed individual 
cells where the two dots had always been on focal planes within the ac- 
quired z-stack during observation. We manually tracked each signal (both 
CFP and GFP) in individual cells and chose the brightest four pixels (2 x 2) 
after the z-stack had been projected to a two-dimensional image at each 
time point. The signal intensities of the selected four pixels along the z-stack 
were summed using Volocity, after background subtraction. The back- 
ground was defined as the most frequently observed signal level among 
individual pixels in images. The intensities of CFP and GFP dots were nor- 
malized in relevant graphs (Fig. 1 , C and D; Fig. SI , B-D) so that their 
intensities were in the range of 80-1 20 before replication and their regres- 
sion curves (see below) did not overlap before or after replication (for clear 
presentation). The nonlinear regression curve following the increase in CFP 
and GFP signal intensities was drawn with Prism 5 (GraphPad Software), 
fitting to curves 



Y = a + , 

where a, b, and LogEC 50 (mid-replication time of each dot) were derived 
constants. The calculated LogEC 50 from the nonlinear regression curve al- 
lowed us to estimate the mid-replication times at a sub-minute level (to the 
first decimal place). The average replication time of CFP and GFP dots 
were defined as 0 min. Cells, in which the difference between CFP and 
GFP mid-replication times was equal to, or less than, 3 min were further 
analyzed; the centroid coordinates of signals from CFP and GFP dots were 
obtained at each time point using Image-Pro Plus and the distances be- 
tween the two dots were measured in three dimensions (Fig. 1 , C and D; 
Fig. SI, B-D). Cells were evaluated for replication at the same factory 
based on the close localization of dots (<350 nm, >2 min) within a —3 to 
+2 min time window (Fig. 2 A, orange rectangle). 

When the distances between the two dots were analyzed before 
(Fig. 2 B) and after replication (Fig. 2 C), the data within the —3 to +2 min 
window were not included; inclusion of such data would have given a bias 
based on the classification itself between the same vs. different factories. 
Moreover, in some cells, close localization continued beyond this evalua- 
tion time window; such data points were also excluded from the distance 
analyses in Fig. 2, B and C, for the same reason. In Fig. 2 C, we aimed to 
analyze the distances between CFP and GFP dots while a replication fork 
actively underwent replication along the right side of the ARS727 replicon in 
the majority of cells. We estimated that such a situation continued, on aver- 
age, until +1 1 min for the following reasons: (1 ) we assumed that, on aver- 
age, a half of the tetO array completed replication at 0 min; then it took 
~4 min to replicate the other half, given that the fork speed is 1 .5 kb/min 
(Sekedat et al., 2010); (2) it took ~9 min for a replication fork to move 
from the right end of the tetO array to the lowest point between the ARS727 
and ARS728 replicons in the replication profile (see Fig. 1 A); and (3) we 
did not include the last 2 min because another fork from ARS728 might 
move in to the relevant region in some cells in the population. 

To define the contours of replication factories in Fig. 3, A and D, 
surface rendering was applied to the GFP-Pol30 (PCNA) signal using 
Imaris software (Bitplane). The objects (factories), which were less than 
15 voxels (~0.002 pm 3 ) in size, were excluded from the analyses to remove 
the background. The total signal intensity of all factories was calculated 
in individual cells (n = 23). The average total signal intensity measured in 
the population was divided by an average number of forks (302) to obtain 
the average signal intensity for a single replication fork. Using this value, we 
estimated the numbers of forks at individual factories (Fig. 3 D). These 
numbers were then binned to even numbers (e.g., 3-5 binned to 4), as sis- 
ter replisomes (therefore sister forks) are associated (Kitamura et al., 2006); 
thus, the numbers of sister replisome pairs were estimated at individual 
factories (Fig. 3 E). In rare cases (3%), the signal intensity of a factory was 
less than the average intensity of a single fork; such a factory was consid- 
ered to contain zero sister replisome pairs (Fig. 3 E). 

The general expression for the Poisson distribution that k events may 
occur is 



P(X = k) = 




where k = {0,1,2,3,...} and X is the mean of distribution. We revised the 
expression considering that there is no factory without a replication fork. 
We calculated the probability for an event of count k to occur under the 
condition that k = 0 was unobservable, as the following: 

P(X = k X > 0) = — = . 

P(X>0) lU(1-e~*) 

The distribution of the number of sister replisome pairs per replication fac- 
tory (Fig. 3 E) was normalized by the sum of the observable counts (for 
k = 1,2,3,...) and then fitted by this zero-truncated Poisson distribution 
P(X = k\X> 0) (Fig. 4 A). The fitted value of X was 1.60. 

Statistical analyses 

Statistical analyses were performed with Prism Software (GraphPad) by 
choosing Fisher's exact test (Fig. 2 A), unpaired f test (Fig. S2 A), c/ii-square 
test (Fig. 3 F), and Pearson correlation (Fig. 4 F), or with R (http://www 
.r-project.org) by choosing Kolmogorov-Smirnov test (Fig. 2, B-D; Fig. S2 B; 
and Fig. S3 Aii). The null hypotheses in these tests (except for Fig. 4 F) were 
that the samples were collected randomly and independently, from the same 
population. For P value in Fig. 4 F, we calculated the chance that random 
sampling from two groups with no correlation would give such correlation 
coefficient. All P values were two-tailed, and the null hypotheses were rea- 
sonably discarded when P values were <0.05. To represent goodness in fit- 
ting, R 2 was calculated with Prism software (GraphPad) as coefficient of 
determination (Fig. 1 , C and D; Fig. SI , B-D; Fig. 4, A and C). 

Mathematical modeling 

(i) "Particles on a string" model. We use an analogy of two particles tethered 
by a string of length d to describe two sister replisome pairs that are apart 
from each other at a chromosomal distance d (Fig. 4 B). Each sister repli- 
some pair is thereby a particle fixed to the end of the string. Each particle 
is considered to be a sphere with diameter e. We assume the string has no 
stiffness, given that the persistence length of yeast chromatin is short (2.5 kb; 
Dekker et al., 2002) relative to the distance between replication origins 
and between marked chromosome loci analyzed here. The particles per- 
form a random walk within a sphere of radius d/2 in three dimensions 
(illustrated in two dimensions, for simplicity, in Fig. 4 B). If both particles 
come within interaction radius (i.e., the distance between their centers is e 
or less), they associate. We fix the coordinate system at the center of mass 
(midpoint) of the two particles. The system can be in two conditions. First, 
when particles are separated; the energy of the system is then E = 0 
(Fig. 4 B, left and middle). Second, when particles are in close proximity 
and become associated; the energy of the system in this state is E = J, 
where J is a binding energy (Fig. 4 B, right). Therefore, J is negative, mean- 
ing that the particles' interaction is attractive. 

Our aim is to estimate the probability of finding the system in each 
condition — particles separated or particles associated — depending on 
the string length between them. The probability that the two particles meet 
and associate with each other when the system is in thermodynamic equi- 
librium is 



P = 2-2 , (1) 

° nB +nB 

a a 5 s 

where n a and n s are the normalized numbers of states in which particles 
are associated and separated, respectively. B a and 8 S represent corre- 
sponding Boltzmann factors (weighing factors). Each Boltzmann factor, 

tw £/,vl , 

depends on temperature, T, and energy, E, of the system, kg is the 
Boltzmann constant. 

The normalized number of states is derived as follows. As the refer- 
ence frame is centered at the midpoint, states corresponding to the parti- 
cles separated by a distance R lie on a spherical shell of radius R/2. 
Particles are considered to associate once the distance between their cen- 
ters is less than s, i.e., they are within a sphere of radius e/2 around the 
origin. The volume of this sphere is 
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We normalize the number states to the total volume 
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The normalized number of states in which the two particles associate is 
then given by 
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The energy of the system at this association state is minimized, thus E a = J, 
with the corresponding Boltzmann factor 



after replication initiation at origins but before fluorescently labeled loci 
are replicated. 

For this, we needed to estimate the diffusion coefficient D of a sister 
replisome pair. We estimated D based on the diffusion coefficient of a 
marked locus around its replication (Fig. 1 A). To measure this, we synchro- 
nized cells with a-factor and then released them into S phase. We started 
acquiring images of the spindle pole body (SPC42-4mCherry) and a fluores- 
cently labeled locus (on ARS728 replicon] 35 min after a-factor release. 
Images were acquired every 7.5 s for 8 min. This time window covered the 
average replication time of this locus (41 min after a-factor release). Because 
the spindle pole body (SPB) is embedded in the nuclear envelope and less 
mobile than a chromosome locus, it was used as a control point for analyz- 
ing the position of the fluorescently labeled locus. We calculated the mean 
squared displacement (MSD) between SPB and the locus as a function of 
time. Diffusion coefficient was obtained from MSD (in four cells) as described 
previously (Marshall et al., 1 997) and we found D s 0.2 um 2 /min. 

To estimate diffusion time-scale, we can compare diffusion size-scale 
with the distance between the two particles. The typical time for one parti- 
cle to diffuse through a distance L is (Sneppen and Zocchi, 2005) 
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The normalized number of states in which particles are not associated is 
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for a small interaction radius (e << d). The energy of the system when the 
particles are apart is E, = 0, and this Boltzmann factor is 8 S = 1 . 
The association probability from equation (1) is then 
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where 
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is a constant with a constant temperature. Equation (4) describes the proba- 
bility of the two particles, separated by a distance <d, being associated in 
an equilibrium system. We can estimate the parameter A from the grouping 
frequencies of replicons in strain #1 and #2. The function P a (d) (shown as 
P(d) in other parts for simplicity) fits the data well (Fig. 2 A), with the best- 
fitting A = 8.7 x 10~ 6 kb~ 3 (R 2 = 0.99, Fig. 4C). The binding energy of two 
sister replisome pairs is J = k B T InjAs 3 ) = —5. 1 k B T = — 1 2.5 kj/mol, for the 
best-fitting A, e = 90 nm and T= 298.2 K. Here we estimate the diameter of 
a single sister replisome pair from the minimum size of a replication factory 
of ~90 nm (Baddeley et al., 2010). We apply the chromatin packaging 
ratio of 1 0 nm/kb, based on the ratio of measured spatial distances over the 
chromosomal distances between two fluorescently labeled chromosomal 
loci, which matches a reported value (Dekker et al., 2002). The calculated 
binding energy of sister replisome pairs (—1 2.5 kj/mol) is in the range of a 
typical weak protein-protein interaction (Baxter et al., 1 998; Rippe, 2007). 
It is also in agreement with the estimated energy for the association of DNA 
polymerases bound on two replication origins (Marenduzzo et al., 2006). 

(ii) Diffusion time-scale. Because replication factories are formed by 
assembly of replisomes that are actively involved in replication (Kitamura 
et al., 2006), it is thought that sister replisome pairs can be associated with 
each other only after replication is initiated at replication origins. Next, we 
addressed whether diffusion of replisome pairs, assumed in the model 
above, is rapid enough for two replisome pairs to meet with each other 



Because we have two particles diffusing, this is equivalent to one particle 
diffusing with 2D. Hence, we double the diffusion coefficient in equation 
(5), D 2 = 2D. For the maximum distance between the origins in strain #2, 
I a 1.3 unn, corresponding to chromosomal distance of 129 kb. With 
D « 0.2 urn 2 / min, the diffusion time-scale is t L ~ 2 min. In strain #2 it takes 
4-6 min after replication initiation at ARS727 and ARS729 until tetO/ 
lacO dots are replicated (i.e., until replication forks reach the middle of 
tetO/lacO arrays, which are 10-1 1 kb in length). We therefore conclude 
that diffusion is rapid enough for two relevant replisome pairs to come 
together before replication of the dots. 

Because the chromosomal distance between the relevant replication 
origins in strain #3 (ARS727 and 731) is larger (1 86 kb) than those in strains 
#1 and #2, diffusion may not be rapid enough for two replisome pairs 
(generated at ARS727 and 73 1 ) to meet before the replication of two fluor- 
escently labeled loci in strain #3. For this reason, the grouping probability 
of the replicon pair in strain #3 was not considered in fitting the parameter 
A (Fig. 4 C). Nonetheless, after A was fitted using the other data points, the 
model was in good agreement with the data point from strain #3 (Fig. 4 C). 

(iii) In silica distribution of sister replisome pairs in replication factories. 
With the "particles on a string" model (Mathematical modeling i), we could 
estimate grouping probabilities of neighboring sister replisome pairs as a 
function of chromosomal distance between them. We then aimed to extend 
our analysis of the replicon assembly from the region on chromosome VII 
to the whole genome context in budding yeast where the positions of ori- 
gins have been well defined by various studies (Nieduszynski et al., 2007). 
For this, we required information about locations of sister replisome pairs 
along all chromosomes, during replication in individual cells. How can we 
obtain such information? 

The replication timing profile of the yeast genome, which was ob- 
tained from population-based studies (Raghuraman et al., 2001; Yabuki 
et al., 2002), helped our understanding of the global temporal organization 
of replication. However, these data provide only the average replication 
timing in the relevant cell population, whereas replication timing is actually 
different from cell to cell. For instance, the same replication origin does not 
fire at the same timing in different cells. Intriguingly, if replication timing is 
obtained in individual cells, the average timing from such dataset could re- 
capitulate replication timing obtained from a population-based analysis 
(Czajkowsky et al., 2008). Importantly, the converse approach is also pos- 
sible. For example, we can set the parameters that determine the behaviors 
of replication origins, such as (1 ) the fraction of cells where an origin was 
competent to fire, (2) mean activation time of an origin, and (3) distribution 
of the origin activation time (a standard deviation, assuming the Gaussian 
distribution). Then, these parameters can be determined so that they reca- 
pitulate the replication profile obtained from a population-based analysis 
(de Moura et al., 201 0). In this approach (de Moura et al., 201 0), it was 
assumed that the speed of replication fork movement along a chromosome 
was reasonably constant (1.5 kb/min; Friedman et al., 1997; Sekedat 
et al., 201 0). Such a mathematical approach was applied first to recapitu- 
late replication timing profile of chromosome VI (de Moura et al., 2010) 
and then further extended to recapitulate replication timing of the whole 
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genome (Retkute et al., 201 2). Once the above parameters were determined, 
we could simulate replication progression in individual cells by assigning 
timing of replication initiation at origins, randomly but still based on the 
relevant parameters. 

In the current study, we used this mathematical approach on two oc- 
casions. First, using this approach, we obtained the distribution of the fork 
numbers on the whole genome in individual cells over the duration of DNA 
replication, as shown in Fig. 3 C. The average number of forks at the peak 
of replication (when the total fork number became maximum) was subse- 
quently used to estimate in vivo distribution of sister replisome pairs in repli- 
cation factories (Fig. 3, D and E). Note that, in Fig. 3 C, we aimed to 
estimate how the number of forks in individual cells changes during S phase. 
For accurate estimation, we needed to consider temporal variation in S phase 
entry among cell population after a-factor treatment and subsequent release. 
We estimated that cells entered S phase with 4-min variation (in standard 
deviation) based on the following result: Globular signals of PCNA repre- 
sented replication factories and its first appearance could define S phase 
onset (Kitamura et al., 2006). We investigated this timing in individual cells 
(n = 56) after a-factor treatment and subsequent release, and found that 
timing of S phase onset varies with 4-min standard deviation. The standard 
deviation was reproducible in repeated experiments. 

Second, we intended to obtain in silico distribution of sister replisome 
pairs in replication factories. For this, we ran the simulation in one million 
cells and took snapshots of replisome positions on chromosomes at the peak 
of replication (Fig. 4, D and E) after cells had entered S phase with 4-min 
variation (see above). Based on these snapshots, we determined whether 
adjacent sister replisome pairs were grouped in the same factory or not, 
based on the chromosomal distance between them and corresponding prob- 
ability of grouping (Fig. 4 C and Mathematical modeling i), as follows. Let us 
designate each sister replisome pair along a chromosome as A, B, C, D, . . . etc. 
in the order. To determine whether A and B are grouped to the same factory, 
we drew a uniformly distributed random number between 0 and 1 (1, but 
not 0, is included), which was then tested against the distance-dependent 
value P(d) of those two pairs. If the random number was below P(d), the 
pairs were assumed to be part of the same factory. Next, we examined the 
association of pairs B and C in the same way. If A and B were in the same 
factory and if B and C were in the same factory, then we concluded that A, 
B, and C are grouped together in the same factory. We performed this pair- 
wise clustering of adjacent sister replisome pairs into factories in the right- 
ward direction along each chromosome. Nonetheless, we confirmed that 
clustering in the left direction gave a very similar result (Fig. S3 C). 

In this study, we assumed that sister replisomes were always associ- 
ated with each other during replication of a relevant replicon. This assump- 
tion was based on our previous result that sister replisomes were associated 
in vivo in most of the cells (Kitamura et al., 2006). We also assumed that, 
when two replisome pairs had encountered upon completion of DNA be- 
tween them, one sister replisome in each pair disappeared, leaving the re- 
maining two replisomes (which originally belonged to two different pairs) 
associated and allowing the new pair to undergo replication. This assump- 
tion is consistent with a low energy state of associated replisome pairs; 
i.e., once two pairs become associated, we can expect that they stay associ- 
ated for a while (see Mathematical modeling i). Nonetheless, in the above 
mathematical simulation, we observed a low number of cases where one 
replisome was present without its sister. This happened when one repli- 
some completed replication at the end of a chromosome (which is linear) 
while its sister was still engaged in replication. This led to generation of a 
small number of replication factories containing odd numbers of replisomes 
(Fig. S3 C). However, for a direct comparison of the distribution of sister 
replisome pairs in factories obtained from in vivo and in silico data, we 
partitioned factories with odd numbers of forks (replisomes) proportionally 
to the nearby categories with even numbers of forks (e.g., factories with 
three forks were recategorized to those with two and four forks proportion- 
ally to their factory numbers). 

In the above mathematical modeling, we assumed that replisome 
pairs A and C could only associate when both A/B and B/C associate. In 
other words, we considered association between immediate neighbors but 
not between others. It was actually not simple to consider direct association 
between A and C because we needed to consider the presence and ab- 
sence of A/C association separately depending on whether A/B and B/C 
association was present or not. In practice, it was impossible to compute all 
possible cases for all possible associations. Nonetheless, our approach is 
justified only when A/C association is relatively low compared with A/B 
and B/C association. We tested this in a simplified case where the chromo- 
somal distances between A and B and between B and C are both d. The 
ratio of the A/C association probability to the A/B and B/C association 
probability was calculated as follows. 



2P(d)' 

When a" = 36 kb (the median chromosomal distance between two neigh- 
boring replisome pairs; see Fig. 4 E) and A = 8.7 x 10~ 6 kb~ 3 , this ratio 
was 0. 1 7. Thus, it was ~6 times more likely to observe A/B and B/C asso- 
ciations than A/C association. 

Online supplemental material 

Fig. SI and S2 provide results supplemental to Figs. 1 and 2, respec- 
tively. Fig. S3 A provides results supplemental to Fig. 3 F and analyses 
the direction in which replicated DNA was reeled out of a replication fac- 
tory. Fig. S3 B provides results supplemental to Fig. 4 C. Fig. S3 C gives 
supplemental information to Mathematical modeling iii in Materials and 
methods. Online supplemental material is available at http://www.jcb 
.org/cgi/content/full/jcb.20 1 306 1 43/DC 1 . 
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